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Abstract 

We analyze stationary accretion of selfgravitating gas onto a compact center within general-relativistic 
radiation iiydrodynamics. Spiierical symmetry and thin gas approximation are assumed. Numerical inves- 
tigation shows that transonic flows exist for small redshifts and they cease to exist for high redshifts and 
high luminosities. There exist two branches of flows (subsonic or supersonic) that originate at a bifurcation 
point and that embrace the set of subsonic solutions. The morphology of the set of subsonic solutions is 
essentially independent of redshifts and flows that belong to their boundary provide estimates of the gas 
abundance of subsonic solutions. It appears that prescribed boundary data guarantee uniqueness only of the 
bifurcation point, and that the latter has maximal luminosity. 
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INTRODUCTION. 



Consider a general-relativistic system — a compact core immersed in a steadily accreting self- 
gravitating gas. The gravitational binding energy of the infalling gas can be converted to a radia- 
tion. Assume that an external distant observer can measure total luminosity, asymptotic tempera- 
ture and redshifts of the radiation. Let be known: the total (asymptotic) mass of the system and the 
physics of the mixture of gas and radiation. Then it would be natural to ask: what mass is within 
the compact body? Alternatively, the mass of the core would be known and the total mass would 
require determination. 

The main goal of this paper is the numerical investigation of this problem for stationary flows. 
We assume spherical symmetry and adopt thin gas approximation in the transport equation. It is 
already known from studies of newtonian radiation hydrodynamics [1] — [3] that supersonic flows 
are generically not fixed by total luminosity, asymptotic temperature and redshift. To each set of 
such data there can correspond two solutions with difl^erent gas abundances. Changing luminosity 
one obtains two curves, on the luminosity-(gas abundance) diagram, that originate at a bifurcation 
point. This point is unique, for given boundary data. General-relativistic supersonic flows with 
small redshifts are similar to newtonian ones in that they also branch from a bifurcation transonic 
flow. In the case of high redshifts supersonic general-relativistic flows can be absent. A similar 
picture appears in transonic flows of perfect gases, newtonian or general-relativistic, without radi- 
ation. In this case boundary data can consist of the mass accretion rate and the asymptotic speed 
of sound [4] and the only unique solution — a branching point — corresponds to the maximal 
accretion. 

Accretion systems with subsonic flows are not determined by the data described hitherto. One 
needs additional information, for instance the asymptotic gas density, in order the specify the so- 
lution completely. We discover, however, an interesting fact valid in the newtonian case and in the 
low-redshift regime of general relativity: transonic flows encompass, on the luminosity-(gas abun- 
dance) diagram, the set filled with subsonic flows. Therefore the two transonic branches provide 
estimates of the mass abundance of corresponding subsonic solutions. In particular, numerical 
analysis suggests that the most luminous flow is supersonic. This picture is valid in the newtonian 
level and also in the general-relativistic case, for small redshifts. If redshifts are large, then the 
boundary of the set of subsonic solutions may consist of transonicor subsonic flows, but it is re- 
markable that the shape of the set of subsonic solutions is only weakly dependent on redshifts. In 
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particular, the flow with maximal luminosity is unique. 

This investigation can be useful in the analysis of two important issues. There is the question of 
identification of the so-called Thorne-Zytkow stars [51], that consist of a hard core and overblown 
atmosphere. They are conjectured to result in the merger of a main sequence star with a neutron 
star or a black hole. If the core consists of a neutron star, then its mass is roughly known. Re- 
sults of this paper show that one can estimate the total mass by measuring luminosity, asymptotic 
temperature and redshifts. Another interesting application would be to distinguish compact stars 
(neutron stars or gravastars i^) from black holes, but the present analysis would require further 
elaboration. Within the scenario investigated here it does not seem feasible to distinguish between 
a black hole or a gravastar, but the investigation of stability can possibly give further information. 

The organization of the rest of this paper is following. Section 2 presents spherically symmetric 
equations of radiation hydrodynamics. The next section explains the concept of quasistationary 
solutions. The interaction of gas and radiation is treated in the thin gas approximation [|7|]. The final 
form of required equations is displayed. Boundary conditions are described in Section 4. The next 
Section brings a discussion of boundary conditions. In particular, we explain the relation between 
the binding energy of collapsing fluids and radiation redshifts. We demonstrate in Section 6 that 
supersonic solutions constitute a one-sided boundary for the set of subsonic solutions in newtonian 
test fluids. Section 7 shows main results of this paper — relations between transonic and subsonic 
flows in the general-relativistic case. The last section contains a summary. 



EQUATIONS OF GENERAL-RELATIVISTIC HYDRODYNAMICS. 

We use comoving coordinates t,r,0 < 6 < n,0 < cp < In: time, coordinate radius and two 
angle variables, respectively. The metric can be chosen in the form 

ds^ = -N^df + adr^ + R^dQ}. (1) 

7? in ([T]) is the area radius. The infall velocity of gas is equal to t/ = and it is related to 
extrinsic curvatures of the Cauchy hypersurface t = const, 

{tJ:K-K;)R = 2U. (2) 

tr^ is the trace of the extrinsic curvature and K'^ is its radial-radial component. 

The energy-momentum tensor reads T^y = T^^ + T^^, where the baryonic part is given by 
^/fv = (P + PWf^Uy + pg^y with the time-like and normalized four- velocity U^, Uf^U^ = -1. The 
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radiation part has only four nonzero components, = -p^ = -T^^ and = T^^, which is 
consistent with the so-called thin gas approximation (see next section). 

A comoving observer would measure local mass densities, material p = T^^^U^Uv and radia- 
tion pe, respectively. The baryonic current is defined as j'^ = PqU^, where po is the baryonic mass 
density. Define as the unit normal to a centered (coordinate) sphere lying in the hypersurface 
t = const and k as the related mean curvature scalar, k = |V,n' = -qfprR- The comoving radiation 
flux density reads j = U^rfNTf I = NT"}^ I Vs. 

We assume a polytropic equation of state for the baryonic matter, p = Kp^^ (K and Y are 
constants) and the relation 

p=Po + h, (3) 

where the internal energy h is easily shown to be equal to p/(r - 1). The equations of motion 
consist, in the spherically symmetric case, of three Einstein equations, of two constraints [|8|] 

(Rk') = -R{Sn (p+Pe) + ^iKf) + + ^(tr i^)^ + tr KK'^, (4) 
ya R R 

and one dynamical equation 

+%nNTl + ^A^(tr Kf - ^ tr KK^. (6) 
The baryonic current is conserved, 

V,/ = 0. (7) 

There are four conservation equations V ^Ty^ = -V^^Ty^ = Tv (here v = 0, r). The quantity Tv is 
called the radiation force density and it describes interaction between baryons and radiation. Its 
radial component will be written as Tr = k^F/, Fr is defined later. The formulation of general- 
relativistic radiation hydrodynamics presented here agrees with that of Park [9] and Miller and 
Rezzola jlOl]. and (on a Schwarzschildean background) with Thorne et. al [11]. One can solve 
formally both constraints, arriving at ||8|] 



K; = dRU-AnRj, 
k 



tiK = ^dK(UR^)-4nR^-. (8) 
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Here m(R) is the quasilocal mass 

m(R) = m{R^) - An j drr^ (p + + -^j . (9) 

The integration in Q extends from R to the outer boundary 7?oo of the accretion system. m(Roo) can 
be equal (or arbitrarily close) to the total asymptotic mass M. The contribution to the mass coming 
from the exterior of Rco can be neglected. Comoving coordinates can be understood as a choice of a 
particular integral-type gauge condition lll2ll . In what follows we will use the comoving spacetime 
foliations, with the time t but often with the areal radius R instead of the comoving radius r. The 
parametrization (t, R) can be interpreted as corresponding to an observer at rest at R. 

One can choose an alternative set of coordinates in a different spacetime foliation, in the so- 
called polar gauge tr^ = (no summation), with the time ts = tsiR, t) and the areal variable R. 
We do not employ these variables here (but see a remark below). 

A simple but lengthy calculation shows that the local mass changes according to the following 
rule 

dtMR) = idt-NUdR)m{R) = An{NkR^{^ii^+(^^'^ + lpEj^ +A^, (10) 

where is the value of -AnNUR^ {p+ Pe + ^) at i? = R^. It is interesting to note that the 
expression AnNkR^ ^7 |l -I- ( j j j -1- 2i7p£ /fcj represents the radiation flux measured by an observer 
located at R in the polar gauge foliation. The mass contained in the annulus (i?,i?oo) changes if 
the fluxes on the right hand side do not cancel. In the case of quasistationary flows the mass is 
approximately constant. 

QUASISTATIONARY FLOWS IN THE THIN GAS APPROXIMATION. 

We will say that the accretion process is quasistationary if all relevant quantities measured in 
the rest frame are approximately constant during time intervals much smaller than certain charac- 
teristic time scale T. In analytical terms, we assume that dt^X = (dt - NUdujX = for X = po, p, 
j, U... 

Under quasistationarity assumption the evolution equation Q can be written as 

^dU k^dN m(R) , , 

From this one easily obtains 

N = kexp^-An J dr^_[p + p + 2pE + (12) 



The local accretion mass rate M is defined as 



M = -AnUR^pQ. (13) 

The baryonic current conservation equation (|7]) takes the form BrM = -I6n^ jR^^\ thus the local 
mass accretion rate M can change, but the quantity 



M = -AnUR^po - (Anf drr' j"^ (14) 

Jr k 

remains constant, = 0. Notice that M(i?oo) = M- 

As the characteristic time scale one can choose the quantity related to the runaway instability, 
T = M/M. One can obtain a rough estimate j^M < CM^ where C is a constant (see dll] for 
the corresponding derivation in the newtonian case). Let Mq be the initial mass. Then M < 
1/ (1 - CMoO- If M/M ^ t then M « Mq and if ? « T then M » Mq; thus the time T sets the 
time scale for the runaway instability. 

The radiation force density has only one nonzero component. It is a simple exercise to show 
that the conservation of M and relation ([3]) imply, for polytropic gases, the vanishing of the zeroth 
component of the radiation force, Tq. Therefore the two related energy and radiation energy 
balance equations read 

NU— +N tr Kip +p) = 
dR 

UR^dR^ ^ ' R^NdR^ ■' ' ■' k 

The two other energy-momentum conservation equations are displayed below. The relativistic 
Euler equation is given by 



„3 .PO 



\ dN ^ ^ dp ^ 



and the transport equation reads 



JLA. (U^R^j) + (n'r^p,) - SnRf- = -NkFr. (17) 

UR^dR^ ■'' R^NdR^ ^ ' ■' k ^ 

There enters an important phenomenological assumption that can be easily expressed in terms of 
quantities related with comoving coordinates. Namely we assume the so-called thin gas approxi- 
mation ilzI] 

Fr = Kpoj- (18) 
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The only direct interaction between baryons and radiation is through elastic Thompson scattering. 
a: is a material constant, depending in particular on the Thompson cross section cr, k = cr/ {4nmpcy 
c is the speed of light and nip is the proton mass. 

In summary, the complete set of equations would consist of eqs ([8]), Q and (fT2l - [r8l) . 



FINAL EQUATIONS. 

It is convenient to express all relevant quantities in terms of the speed of sound a, given by = 
^, because we will search for flows possessing transonic points. We find that it is computationally 
expedient to replace the radiation energy balance equation (the second of eqs (fTSl) ) by the total 



MN— — ^— + 2MN— = 4nR^jNk(l + I + C. (19) 



energy conservation 

F - 1 - Po \ 

Equation (fT9l ) is a direct consequence of quasistationarity and the equation (flOl) . The constant C 
is the asymptotic flux (i.e., flowing through the sphere of a radius Rco) in (fTOl) . 
Furthermore, write down eq. (fTTI) as 

— (R^j) = -KkNjpo T — [R Pe] - 

dU dN f 

INj— - 2kpE— + SnNR^-. (20) 
dR dR k 

and replace the term {R^j) second of equations (fTSl) . After some algebra one arrives at 

' 2m{R)\ N d / , ^ J dU 

2k(JU-kpE)^ + inNR{j^-jPEjy (21) 



Below we display a number of relations, that can be easily obtained from the equation of state and 
([3]), namely 

p_ _ F - 1 

Po F F - 1 - ' 
p+p ^ F-1 
Po F - 1 - ' 



p+p 



5Rln|^^j = -dR\n{r-\-a^), 



a V-' r r-1 



Po = Pooo — r- (22) 



It is useful to insert (|22l) into Eqs. ([14]— [TS]). One obtains following equations: 



) the gas energy density conservation equation (the first of eqs. (fT5] )) 



cl 7 r - 1 - 
— In a = 7— X 

dR a^-f 

f)) 



^iii-fnin (23) 



ii) the baryonic mass conservation 

dU U d ^ , 2U AnRj 

— = T — Ina^ + (24) 

dR Y-\-a^dR R k ^ ^ 

iii) the equation for the lapse 

— =N\Kj + — In(r-l-a^) . (25) 

dR [■' r-l dR ^ >j ^ ^ 

Eq. (ESI) follows from ([M]) and (ED). 

Equations ([191), dlD and (|23l— [251) constitute, with k and m{R) given by ([8]) and ((9|), the complete 
model used in numerical calculations reported in next sections. A remark is in order. It is clear 
from the inspection of ([231) that if = (the speed of sound equals the spatial length of the 
infall velocity) then the expression ^ (f - 2t/^ + AnR^ i^pE + p + jjJ^-Kj(^l - must vanish. 
There are four different ways of passing through the transonic point (similarly as in the newtonian 



analysis in ill BID and only one of those corresponds to the accretion. 



BOUNDARY CONDITIONS. 



The overall picture of the system is as follows. A ball of gas is enclosed by a sphere Soo of 
a radius i?co and connected, via a narrow transient zone filled with baryonic matter and radiation, 
to the Schwarzschild vacuum spacetime. It is clear that by careful arrangement of data the mass 
within the transient zone can be negligible. Therefore we assume that the asymptotic mass M is 
equal to m (7?oo) (see Eq. 

Boundary conditions at the outer sphere 5co are needed for the radiation quantities j, pe, the 
mass accretion rate M, the square of the speed of sound and the baryonic mass poo- We assume 
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that al, » M/Roo » U^; the second inequality means that infall velocity is much smaller than the 
escape velocity. These inequalities guarantee the fulfillment of the Jeans criterion for the stability 



(see a discussion in 



and studies of stability of accreting flows in newtonian hydrodynamics 
Ill7l1 ). suggesting the stability of solutions. 

One can derive, after some algebra involving manipulations of equations and (I2TI) . the 
approximate equation j^{(Pe - j) R^) ~ in the asymptotic region. This means, taking into 
account the fact thati?oo can be arbitrarily large, that one can safely assume jco = Pe„ ■ Furthermore, 
the total luminosity is with good accuracy given by Lq = AnR]^jea and it must be related to the 
accretion rate by the formula 

Lq = aM^. (26) 
The coefficient a determines the relative binding energy. We assume 



(27) 

kiRo) ^ Ro 

where Rq is the outer radius of the hard core of the system. The last two formulae can be justified 
by two arguments. First, in the nonrelativistic limit one gets a = \(f> (Rq) |; a is equal to the absolute 
value of the newtonian potential on the surface of the hard body. It is clear now that (|26|) is 
just the statement that all available binding energy is transformed into radiation, and that there 
is an implicit assumption that the heat capacity of the core is negligible. Second, the condition 
of stationarity implies the existence of the approximate time-like Killing vector. By employing 



standard reasoning nlm . one arrives at the two formulae (|26l ) and (1271 ). Thus a can be regarded 
as a proper binding energy. Let us remark that a gives the standard measure of the gravitational 
red- or blue-shift. If stationary observers detect coq at Rq and co at infinity, and l/co 2M (the 



geometric optics condition — see 111911 for a discussion) then lo = (\ - a) coq. 

In conclusion, boundary data consist of the binding energy coefficient a, total luminosity Lq 
and the asymptotic speed of sound a^. Not only Lq but also the two remaining quantities can be 
in principle determined from observations: a from the measurement of the highest redshift and 
al, from the asymptotic temperature. Then jao = Pe„ = Lq/ (AttRI,^, and the mass accretion rate 
M = Lola. These data in fact specify transonic ffows up to, possibly, a bifurcation; for given data 
there can exist two solutions. In the case of subsonic flows another boundary condition is needed, 
for instance the baryonic mass density pco- We show later that transonic flows can give bounds 
onto some characteristics of relevant subsonic solutions. 
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SUBCRITICAL VERSUS CRITICAL: LESSONS FROM NEWTONIAN HYDRODYNAMICS OF 
TEST FLUIDS. 



In the Bondi model [11311 the selfgravity of gases is neglected. The relevant equations can be 
obtained from these presented above in the following way: assume the test gas approximation, 
a = j = pE = 0, k = N= l and r-l-a^«r-l. In this approximation p = pq. The whole 



problem reduces to two algebraic equations 



M = C (28) 



— + = (29) 

2 r- 1 r- 1 ^ ^ 



From (|28l) one has U = Since now p = poo y-^j ' , one obtains 



t2 ^fi 



2 ^Th 



U^ = 2f,[^j , (30) 
where /3 = -^^TpT- Insertion of dSO]) into ^ yields 

R^\a^l r-1 R r-1" ^ ^ 

Thus newtonian transonic flows have this interesting property that the quantity /3 is given by a 
simple analytic formula involving boundary data: 



I 2 



5-3r 



that means that, given M and a^o, the density poo and p is also specified. One finds from Eq. (|3T1) 
that 

da^ U^a^ 1 



and there exists (from the implicit function theorem [|14l1 ) a local solution a^(J3). One infers from 
(l33l) that outside the supersonic sphere the speed of sound decreases with the increase of /3. Taking 
that into account, one concludes from the inspection of (l30l) . that the infall velocity increases with 
the increase of /3. Therefore subsonic flows, for which < a^, can exist only for values of the 
parameter /3 smaller than yS^- That means that, given a transonic and subsonic flows with the same 
data Geo, M, the asymptotic mass density of the subsonic flow must be bigger than that (pcco) of the 
(test fluid) transonic flow, 

(34) 
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Notice that the asymptotic mass density of transonic flows can be represented as pcoo = fic/M; 
therefore the mass density is now completely specified by the other data. Another way of rephras- 
ing (|34l ) is to say that test fluid flows are more eflicient than subsonic ones in the sense that a given 
mass accretion rate demands less gas in the first case than in the other. See also another derivation 



in u5\\ . One can infer from this description that the branch of test fluid transonic flows embraces 
from below the set of subsonic solutions, in the (mass accretion rate)-(gas abundance) diagram. 
The same conclusion holds true for the Shakura model [16], assuming the test fluid approxima- 



tion. Indeed, in this case one has instead of Eq. (|32l) the following yj] 

£_ l^y-^ _M_ _ai_ 

R^\a^l r-1 R r-1' 

where only and /3 depend on the mass parameter poo. Here (f>(Ro) is the newtonian potential 
on the surface of the compact core and G denotes the gravitational constant. One can find the 
dependence of as the function of poo in a similar way as before. Thus a supersonic flow with a 
given luminosity can have less gas than a subsonic flow of the same luminosity. 

In particular, again we observe that subsonic flows lie above the supersonic (test fluid) branch, 
in the luminosity-(gas abundance) diagram. Numerical studies show more, that the set of subsonic 
solutions has a parabola-shaped boundary that consists only of transonic solutions. The edge of 
the parabola has maximal luminosity. 

NUMERICAL RESULTS. 

Below we will study how much of the newtonian picture drawn in the preceding section is valid 
in the general-relativistic case. 

The equations are put in the evolution form, see eqs. (fT9l) . (|2T)) . (|23H25I) . We start from the outer 
boundary Rco and evolve inwards until the equality 



NiR) 2m (R) 

is met. The corresponding value of the radius is denoted as Rq and it is regarded as the size of the 
inner core. It appears that the 4th order Runge-Kutta method fails almost immediately and for that 



reason we employed the 8th order Runge-Kutta method pOl] 
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We assume required boundary data joo = P£„, the mass accretion rate Moo, the speed of sound 
a^, the parameter a and the baryonic mass density poo- Inequalities Ul, M/Roo are 
obeyed, since this ensures the Jeans length of the configuration to be much larger than R^o, which 
in turn suggests stability of the configuration (see a discussion in [1]). Solutions are obtained by 
the method of shooting. Given a^, a and Lq = aMoo, one should vary the asymptotic baryonic 
mass density pooo. We find that solutions exist for a finite range of bi < pooo < ^2- The change 
of the luminosity Lq, while keeping constant and a, results in another segment of solutions 
with baryonic densities b[ < pooo < Solutions corresponding to the extremal baryonic densities 
bi(Lo), b2{Lo) are shown in Figs 1-3 as the bifurcation curves. Each point on the first three figures 
(1-3) within the region enclosed by bifurcation curves represents a solution with a specific value 
of the asymptotic baryonic density and fixed boundary data (a^, a and Lq = aMoo)- 

The numerical integration is straightforward with the important exception of transonic solu- 
tions. It is clear from Eq. (I23l) that the flow becomes critical at a sonic point and if = U^/k^ 
then 

The numerical strategy for finding transonic flows is as follows. For a density pooo chosen at 
random one either obtains no solution at all or a subsonic solution. Using the bisection method 
one can obtain a boundary of the solution set, later on called the bifurcation curve. This search 
process can be automated and it works well for all values of the parameter a. For small parameters 
a one gets convincing numerical evidence that the bifurcation curve consists solely of transonic 
solutions. 

When a is close to one, then the automated search produces a bifurcation curve, but the question 
whether it consists of transonic flows has to be studied in a more detailed way. If a is significant 
then the numerical problem becomes quite sensitive on tiny deviations - of the order of lO^^Vo - 
from the right values of the asymptotic mass density. The bifurcation curve is found within some 
margin error and that error would in many cases be larger than lO^^Vo; thence automated search 
becomes inconclusive. Investigating collected data one can in some cases determine the character 
of a solution that lies on the bifurcation curve. 

Another subtle problem is evolving Eq. (l23l) in the vicinity of the sonic point. The analytic 
reason is due to the fact that the sonic point is a critical spatial point at which coalesce four diff'erent 
solutions: two accretion branches and two wind branches, in each case one inside and one outside 
of the sonic sphere. This agrees with the well known feature of the standard Bondi accretion of test 
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fluids [|l3h . The accretion flow solution consists of two branches, that existing outside of the sonic 
sphere and the other that bifurcates inward from the sonic point. The accretion branch is unstable 
beneath the sonic sphere. From the numerical point of view when the denominator - jk^ of 
Eq. (|23l) is small (smaller than 10"'^) then the whole fraction is calculated with a large error. For 
that reason there must exist a procedure for checking the value of - U^/k^ at a point and if it 
becomes too small then a regularization method must be implemented. It appears that it is enough 
to do this regularization only during one step while passing through the sonic point. 

We assumed following asymptotic parameters. We choose Mq/M = 5.95496 x 10"^, where 
Mq is the Solar mass. Let us remark here that, while a definite choice is needed by the nature 
of numerical calculation, there is nothing pecular in the above data. One can repeat this analysis 
assuming that the total mass is of the order of the Solar mass. 

The parameter k = cr/ (4nmpcj, as we explained earlier. In this paper we adopt the traditional 
choice of units G = c = 1 and supplement it by the scaling M = I. This leads to the value of the 
K = 2.1326762 x lO^^ (Mq/M), that is k = 1.27 x 10^^ The size of the system is Roo = 10^ The 
speed of sound is given by = 4 x lO"'* in all numerical calculations that are described below. 
The Eddington luminosity reads Le = AttM/k = 9.9847 x 10"'^ 

Figures [BS show accreting solutions on the diagram luminosity-(mass of the central core). 
Each point within the set embraced by two curves (bifurcation curves) corresponds to an accretion 
solution. Solutions are absent outside of this region. Typically, for small binding energies (that 
is, small a) and small luminosity Lq, there exist two accreting solutions possessing sonic points, 
with asymptotic densities pocoi and pocoi- There appear subsonic flows for each poco 6 (pooouPoooi)- 
Thus the two bifurcation branches of transonic flows embrace a set of subsonic flows. Accreting 
stationary flows are absent above the bifurcation point, which maximizes the luminosity. All that 
is illustrated in Figured] 

Two forthcoming figures demonstrate that with increase of the parameter a the shape of the 
set of subsonic solutions does not change significantly. Its boundary, however, can consist both 
of subsonic or transonic flows. We call corresponding solutions lying on the bifurcation curve 
as extreme. Figure |2] is done for a = 0.5. Figure |3] presents the bifurcation curve for a = 0.9. 
The bifurcation point with the luminosity Lq = 0.31 130L£ on the extreme curve corresponds to a 
subsonic solution (see Fig.[8]for the behaviour of its speed of sound and of the infall velocity). A 
transonic solution exists for Lq = O.IL^, a = 0.9 on the lighter branch, but the partner lying on the 
more massive branch is a subsonic flow, as shown in Fig. |4l 



13 




FIG. 1: Bifurcation curve for small binding energy, a = 0.0025. Two branches of transonic flows encompass 
the set of subsonic flows. The abscissa shows the luminosity and the ordinate shows the mass of the compact 
core. 

Figure |5] shows how squares of the speed of sound and of the spatial velocity U^/k^ depend 
on R for different gas abundances. The flow with greater gas abundance (hence with a relatively 
lighter compact center) possesses a sonic point that is closer to the center than in the other case. 
Intuitive explanation is that for greater gas density the radiation pressure is bigger and prohibits 
quick falloff; the infall velocity can approach the speed of sound only close to the gravity center. 

Figure [6] sketches the behaviour of the speed of sound and U^/k^ for one of the subsonic solu- 
tions, with identical boundary data - but different baryonic densities - as for the flows depicted in 
Fig. [51 

The next figure shows characteristics of the bifurcation solution for the binding energy param- 
eter a = 0.0025. The striking feature is that the content of gas abundance approaches 0.31. We 
observed earlier that in the pure hydrodynamic accretion ^ the maximum accretion rate occurs 
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FIG. 2: Bifurcation curve for intermediate binding energy, a = 0.5. Two bifurcation branches encompass 
the set of subsonic flows. The abscissa shows the luminosity and the ordinate shows the mass of the compact 
core. 

n 

when the gas abundance is equal to 1/3. That is valid both in the general-relativistic [4] and newto- 
nian [21] case. In the Shakura model [1] one can analytically prove that at the bifurcation point the 
luminosity is maximal and the gas abundance must be smaller that 1/3. The inspection of Figures 
1-3 shows that this property is satisfied also in the general-relativistic case. Here brightest flows 
coincide with bifurcation points. Their gas abundance is always smaller than 1/3 and decreases 
with the increase of Lq/Le- The last figure displays the dependence of and U^/k^ on the area 
radius R in the case of the most luminous solutions corresponding to a = 0.9. It is clear that it is a 
subsonic flow. 
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FIG. 3: Bifurcation curve for high binding energy, a - 0.9. Two bifurcation branches encompass the set of 
subsonic flows. Abscissa shows the luminosity and the ordinate shows the mass of the compact core. 

CONCLUDING REMARKS. 

The present analysis is fully general-relativistic, with the backreaction effects included. The 
picture emerging here is different from the former investigation of spherical accretion in which the 



space-time geometry has been fixed and therefore backreaction has been ignored [|22ll . The main 
new feature is the existence of an upper limit for the asymptotic baryonic mass density and of a 
massive bifurcation branch. The new striking element is the fact that the brightest configuration 
is unique for given redshift a and asymptotic speed of sound. In the case of low a (which is 
associated with low luminosity) the brightest object flow is rich in gas — about 1/3 of its mass is 
in the gaseous accreting matter 

On the other hand general-relativistic radiating systems with accreting gas behave in a quali- 
tatively similar way to self gravitating newtonian ones for small redshifts (i.e., for small binding 
energies). In both cases one observes the following feature: two arms of the bifurcation curve of 
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FIG. 4: High binding energy, a = 0.9. A pair of extreme subsonic and transonic flows corresponding to the 
same boundary data Lq = O.ILe and = 0.0004 and diff"erent baryonic densities. Abscissa shows squares 
of the speed of sound and infall velocity and the ordinate shows the radius R. A smaller picture inlet in the 
lower figure shows the vicinity of the sonic point. 

transonic flows embrace, in the diagram Lq - m(i?o) (or the gas abundance versus asymptotic lumi- 
nosity Lo), a set of subsonic solutions (see Fig.[T])- For given asymptotic data (M — the total mass, 
Lq — asymptotic luminosity and — the asymptotic speed of sound) there exist two transonic 
solutions and an infinite number of subsonic solutions with intermediate baryonic densities. The 
mass of the central core in subsonic flows is comprised between two bounds of the two limiting 
transonicflows. The bifurcation point, where the two supersonic branches cross and the luminosity 
is maximal for a given (binding energy) parameter a, is unique. 

Some features of this picture change in the case of larger redshifts. This happens at some 
value of the parameter a larger than 10"^. The massive part of the bifurcation curve is replaced 
by a curve of subsonic solution (inspect Figs |2] and |3]). Transonic flows survive only on that 
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FIG. 5: a - 2.5x 10 ^, Lq/Le - 0.00318008. Behaviour of and U^/k^ in two supersonic flows generated 
by the same boundary data but different baryonic denisties. 

section of the test fluid bifurcation branch where most of the mass is comprised in the compact 
core and the luminosity is relatively low. One observes (for instance for a = 0.5 and a = 0.9) 
that as one increases the total luminosity Lq, the transonic flows completely cease to exist on the 
bifurcation curve and they are replaced by subsonic solutions. Nevertheless the set of all flows has 
a similar shape that in the case of small or, as exemplified in Fig.[T]-[3l Subsonic solutions are not 
specified uniquely for given boundary data, as we point out earlier, but the length of the interval 
of allowed values of the asymptotic baryonic density pooo becomes shorter with the increase of Lq. 
The solution corresponding to the maximal luminosity is unique. In particular flows corresponding 
to highest possible luminosities (that can be close to the Eddington luminosity Le if a is close to 
1) are uniquely determined. 

Finally, it is interesting that quasi-stationary solutions of the model considered in this paper 
can have a significant abundance of the gas. Accreting systems with maximal luminosities (in 
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FIG. 6: a = 2.5 x 10"^ LqILe = 0.00318008. Behaviour of and U^jk^ in a subsonic flow generated by 
the same boundary data as in Fig. [5] 

particular close to the Eddington luminosity) can possess even 33% of gas for small redshifts and 
still almost 10% of gas for a = 0.9. It is an open and important question whether this picture is 
valid for generic nonspherical flows. 
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